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Theory  of  Space  Charge  Between 
Parallel  Plane  Electrodes 


I.  INTRODUCTION 

The  theory  of  space  charge  between  parallel  plane  electrodes  has  been  treated 

by  numerous  authors  beginning  with  Child's  work  in  1911 1  and  culminating  in  a 

2 

fairly  detailed  paper  in  1938  by  Fay,  Samuel,  and  Shockley.  The  present  report 
differs  from  these  earlier  publications  in  a  major  respect;  namely  that  the  earlier 
papers  treat  an  initial  value  problem,  while  the  current  work  treats  a  boundary 
value  problem.  From  a  practical  standpoint,  if  similar  ranges  of  parameters  are 
covered,  this  difference  is  of  small  consequence,  since  for  equivalent  families  of 
solutions  the  information  contained  is  equivalent.  F rom  the  mathematical  point  of 
view,  however,  the  treatment  as  a  boundary  value  problem  differs  substantially 
from  the  initial  value  treatment 

The  planar  space  charge  problem  is  one  of  the  simplest  naturally  occurring 
problems  in  electrostatics.  Since  most  problems  in  electrostatics  are  formulated 
by  necessity  as  boundary  value  problems,  it  is  of  interest  to  solve  the  planar  space 
charge  problem  in  this  manner.  The  various  types  of  solutions  are  also  more  easily 
organized  and  classified  with  this  treatment.  Finally,  such  a  solution  serves  as  a 
rare  prototype  which  can  be  solved  analytically  to  solutions  of  more  difficult  problems. 


(Received  for  publication  6  Decembt  ■  1982) 

1.  Child,  C.D.  (1911)  Phys.  Rev.  32:492. 

■  ~  ■— 

2.  Fay,  C.E. ,  Samuel,  A.  L. ,  and  Shockley,  W.  (1938)  Bell  Systems  Technical 

Journal,  17:49. 


2.  FORMAL  CLASSIFICATION  OF  THE  SOLUTIONS  OF  THE 
SPACE  CHARGE  EQUATION 


The  space  charge  equation  for  a  plane,  written  in  MKS  units  is 
2 

e  i-JC  =ne=J(2e/M(v  +  v) \  1 

°d?  1  °  > 


where 


v  =  electric  potential, 

n  =  charged  particle  number  density, 

e  =  electronic  charge, 

J  =  current  density, 

M  =  particle  mass, 

x  =  spatial  coordinate, 

t  =  permittivity  of  free  space, 

vq  =  initial  energy  per  unit  charge  of  emitted  particle. 


Let 


y  =  x/L,  A  =  9/4 


J  V 


■2S 


e  v 
o  o 


TfT  * 


<0  =  V  /  VQ 


where  L  =  distance  between  electrodes. 

Then 

d2«p  =  4/9  A  _ 

dy^  Vl  + 


Integrating  twice,  we  have 


where  C  and  C*  are  constants  of  integration.  Let  (5=0  when  y  =  0.  Then 


i'U?.A_yj..tJgrI-C  (2C  +  16/9  A  Vl  +  <P>1/2  =  ±  y  +  %  (2 C  +  16/9  A)1 12 


16/27  A 


16/27  A 


(5) 


The  common  space  charge  limited  solution  sets  gy  =  0  when  y  =  0.  Here,  we  treat 
the  general  case  by  letting  <p  =  0q  when  y  =  1.  There  are  then  three  distinct  classes 
of  solutions  which  must  be  distinguished: 

(I)  A  potential  minimum  for  0  5  y  -  does  not  exist. 

(ID  A  potential  minimum  at  0  -  yM  s  1  does  exist,  the  corresponding 
value  of  potential  <5^  is  greater  than  -1. 

(Ill)  A  potential  minimum  at  0  5  y^j  S  1  exists;  the  corresponding 
potential  0^  =  -1;  and  the  current  divides  at  y  =  y^. 

We  examine  each  of  these  in  turn. 

For  class  I,  a  single  solution,  with  either  the  +  (0Q  >  0)  or  -  (iPQ  <  0)  sign  is 
continued  from  y  =  0  to  y  =  1.  The  imposition  of  0  =  <0Q  when  y  =  1  then  leads  to 
the  following  equation  for  the  determination  of  C 


4/9  A 


16/27  A 


(2  C  +  16/9  A 


/l  +  0jl/2  =  ±  1  +  -4 1 9  (2  C  +  16/9  A)112 


16/27  A 


(6) 


If  this  equation  is  cleared  of  fractional  powers  it  results  in  a  cubic  in  C, 

2  C3  -  (1  -  B2)2  C2  -  16/27  A  (2  A  +  (1  -  B2)  (B3  -  1  -  A)]  C 

+  <8/27)2  A2  (4  A  -  (B3  -  1  -  A)2]  =  0  (7) 


where 

B  =  V^o 

Eq.  (7)  contains  implicitly  roots  corresponding  to  the  ±  sign  appearing  in 
Eq.  (6);  it  may  also  contain  extraneous  roots  introduced  by  the  clearing  of  fractional 
exponents.  Thus  any  (branch  of)  solutions  for  C  of  Eq.  (7)  must  be  tested  in 
Eq.  (6)  to  insure  that  they  actually  represent  solutions  of  the  original  boundary  value 
problem.  There  are  thus  both  advantages  and  disadvantages  in  working  with  Eq.  (7) 
rather  than  directly  with  the  original  Eq.  (6).  The  advantage  is  that  Eq.  (7)  is  of 
a  well  known  explicitly  solvable  form,  so  that  the  roots  can  be  systematically  found 
and  classified;  the  disadvantage,  noted  above,  that  once  a  root  is  found  it  must  be 
checked  further  to  see  if  it  also  satisfies  Eq.  (6). 


The  general  explicit  solutions  of  Eq.  (7)  are  too  lengthy  to  be  of  value  and  we 
do  not  write  them  down  here;  we  discuss  rather  their  general  properties. 

The  class  I  solutions  may  be  further  divided  as  follows; 

(la)  <Po  2  0  . 

The  +  sign  in  Eqs.  (3)  through  (6)  Is  chosen  and 
^  )  =  16/9  A  +  2C  2  0  . 

dy/y=o 

For  any  given  id^  there  is  a  maximum  value  of  A,  corresponding  to  C  =  -8/ 9  A,  and 
obtained  by  putting  this  value  of  C  into  Eq.  (6) 

A+m=  (B+2)2  (B-l)  .  (8) 

(lb)  (£>o  <  0  . 


The  -  sign  in  Eqs.  (3)  through  (6)  is  chosen  and 


^jy=l  = 


16/9  A  B  +  2C  *5  0  . 


For  any  given  <£>q  there  is  again  a  maximum  value  of  A,  here  corresponding  to 
C  =  -  8/9  A  B  and  obtained  by  putting  this  value  of  C  into  Eq.  (6) 


A^=(l+  2B)  (1  -B). 


(9) 


The  maximum  value  of  A^  is  obtained  for  =  -3/4  (B=  1/2)  and  is  given  by 


(*m) 


=  2. 


MAX 


Class  II,  - 1  -  <0  s  00 
o 


For  this  class,  the  curve  of  (0  vs  y  is  split  into  two  branches  and  Eq.  (6)  cannot 
be  used  directly  to  fix  C.  The  separate  boundary  conditions  for  the  two  branches 
are: 
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Left  hand  branch 


Might  hand  branch 


(-  sign  applies 
y  =0  0=0 

y  =  'VM  °=UAI*ay]M 


y  =  y 


d°\  =  0 


V  =  1 


(4  sign  applies) 

M  0=  %  37  )  =  ° 
”  M 


o =  o 
o 


For  the  left  hand  branch  we  use  Eqs.  (3)  and  (5)  with  the  -  sign  to  obtain 


:c  4  16/9  A  V 1  +  e>Aj  =  0  with  -1  i©  t  0 


( 10a) 


0  =  '  9  (2C+  16/9  A),/2.  0  y  -  1  . 

16/27  A“  M 


(10b) 


For  the  right  hand  branch  we  use  Eq.  (4)  with  the  +  sign  and  obtain 


—  9  AB  'S  (2  C  4  16/9AB)1/2  =  1  -  y. 
16/27 


M 


(11a) 


—  9 --B  -'-r,-  (2  C  4  16/9  AB) 1/2  =  1  .  _4  / 9  A  -  C  (2  C  4  16/9  A) 1/2  .  (lib) 


16/27  A 


16/27  A 


The  potentials  are  given  by 


—  9 -A  +  (2 C  4  16/9  A  Vl  +  o)i/2  =  ±  (y.,  -  y) 

1  C  /  O'?  \  ^ 


16/27  A 


(12) 


where  the  4  sign  goes  with  the  left  hand  branch  and  the  -  sign  with  the  right  hand 
branch. 

The  only  difference  between  Eq.  (lib)  and  Eq.  (6)  is  in  the  signs  on  the  right 
hand  side.  When,  however,  Eq.  (lib)  is  cleared  of  fractional  powers  to  form  a 
cubic  in  C,  these  differences  are  eliminated,  and  once  again  we  obtain  Eq.  (7). 
Again,  clearing  Eq.  (lib)  of  fractional  powers  may  introduce  spurious  roots,  so 
that  solutions  of  Eq.  (7)  must  be  examined  carefully  to  see  if  they  satisfy  both 
Eqs.  (10a)  and  (lib),  and  are  thus  proper  solutions  of  the  original  boundary  value 
problem. 

The  lower  limits  of  A  for  class  II  solutions  are  the  upper  limits  for  class  I, 
namely  A^.  If  all  solutions  were  unique,  the  upper  limits  would  be  found  from  Eqs. 
(10a)  and  (10b),  the  former  with ©M  =  -1  leading  to  C  =  0,  the  latter  yielding 
7M.1/VA.  Putting  these  values  in  Eq.  (11a) 


[a../4 *■]’• 


In  fact,  however,  not  all  pairs  of  values  (A,  B)  do  lead  to  unique  solutions,  and 

Ay  does  not  always  give  the  upper  limit  for  class  II  solutions.  It  does,  however, 

play  a  significant  role  in  the  classification  of  solutions. 

For  some  values  of  we  find  A.+  £  A*  A*  while  for  others  this  is  not  the 
o  L  M  M 

case.  The  transition  occurs  for 


(By  +  2)2  (By  -  1)  =  (By/2  4  l^ 


and  for 


(2  By  +  l)2  (1  -  By)  =  (By/2  4  l^ 


for  <P  >0 
o 


for  <P  <  0  . 
o 


Eqs.  (14)  and  (15)  when  cleared  of  fractional  powers  are  quartics  with  one  root, 
which  clearly  does  not  satisfy  Eqs.  (14)  and  (15)  at  By  =  By  =  1.  If  these  quartics 
are  divided  by  By  -  1  and  By  -  1  respectively,  the  resulting  cubics  are 

9  B2  4  5  B2  -  25  By  -  25  =  0  (16) 

25  B2  4  25  B2  -  5  By  -  9  =  0  (17) 

from  which  we  see  directly  that  By  =  1/By. 

Solving  Eq.  (16)  we  have  for  the  one  real,  positive  root.  By  =  1.  81612  and 
hence  By  =  0.  55062. 

The  corresponding  values  of  are 


=  2.  2983  and  to  |  =-0.69682. 

°/L 


The  maximum  value  of 


Alil  (B  4  2)2  (B  -  1) 


(B  '  4  1) 


is  4/3  and  occurs  for  B  =  4  or  <f>  =15.  t  he  maximum  value  of 

o 


12 


-15/ 16. 


A, 


‘M  (1  +  2  B)~  (1  -  HI 
"T^T 


IB' 


+  1)' 


is  also  4/3  and  occurs  for  B  =  1  4  or  C 

o 

Tlu>  significance  ol  the  separation  of  values  of  o  into  those  for  which 

V  i  O  \  and  O  \  <D  \  will  be  clarified  in  what  follows, 
o  o Jv  o  o  )  i 


3.  DISCUSSION  OF  SOLUTIONS  FOR  THE  INTEGRATION 
CONSTANT  FOR  CLASSES  I  AND  II 


We  commence  our  examination  of  the  functions  C(A,  B),  implicitly  defined  by 
Eq.  c7 ),  bv  an  application  of  Liescartes'  rule.  We  note  that  the  coefficients  of  the 
first  and  second  terms  are  always  positive  and  negative  respectively.  The  coeffi¬ 
cient  of  the  third  term  is  positive  for  small  A  and  changes  sign  when 


AK  = 


<B2  -  1)  (133  -  1) 


while  the  coefficient  of  the  last  term  is  negative  for  small  A  and  changes  sign  when 


±  'M2  2 

A*  =  (B  '  ±  1) 


.Now, 


<B:i/2  +  o-Mjr-iMB3-,) 


B“  +  1 


for  all 

0  <  B  <» 

while 

{nM2  _u2  .  ( n2  -  1)  (H:i  -  1) 

1 

the  equality  holding  only  for  B=  1.  We  thus  have  the  following  sequence  of  varia¬ 
tions  in  signs  of  coefficients  as  A  goes  from  0  to  . 


13 


A 


Signs  of  Coefficients 


+  -  + 

+  -  +  + 

+  -  -  + 

+ 

The  change  in  sign  from  the  second  to  the  third  row  does  not  change  the  number  of 
variations  in  sign  and  hence  we  have  for  the  number  of  ±  roots,  N±  according  to 
DesCartes1  rule, 

N  N 

3  or  1  0 

2  or  0  1 

1  2  or  0 

The  ambiguities  in  N+  and  N  are  resolved  by  examination  of  the  discriminant  of 
Eq.  (7). 

If  we  use  the  common  notation  for  the  coefficients  of  cubics 

C3  +  pC2+  q  C  +  r  =  0  (18) 

2,o 

a  =  q  -  p  /  3 

b  =  2/27  p3  -  1/3  pq  +  r 
then  the  discriminant  is  defined  by 
R  =  b2/4  +  a3 / 27  . 


0  <  A  <  Al 

AL  <  A  <  Al 
A^  <  A  <  00 


0  <  A  <  Aj^ 

al  <  A  <  Ak 
ak  "  A  "  a+l 

At  <  A  <  * 


The  discriminant  of  Eq.  (7)  is  of  eighth  order  in  A,  and  is  factorable  into  the 
following  form: 


r  =  iii>  A3  (A  -  (1  +  J3)3]  |  A  -  (I  -  ID2  ( B  +  b]2  I A  -  (1  -  U)2(H/2  +  l))2 
(27  C 


(19) 


3  3 

Thus,  for  A  <  (1+B)  there  are  three  real  roots  and  for  A  >  (1  +  B)'  only  one.  Let 

A_  =  (1+  B)3;  then  the  numbers  of  roots  are 
K 


N+  N 


0  <  A  <  A’  3 

aL<a<aL  2 

a+l<a<ar  1 


0 

1 


2 


A  <  A  <  » 


1 


0 


With  this  table  and  the  preceding  material  in  hand  we  are  in  a  position  to 
organize,  and  in  what  follows,  to  assign  types  of  solutions  in  the  parameter  space. 
Plotting  the  four  curves  A^j,  A A*  ,  and  in  the  A,  B  (or  A,  tPQ)  plane,  we  find 
the  quadrant  A,Ht  0  divided  into  seven  regions.  Each  region  will  be  characterized 
by  the  type(s)  and  numbers  of  solutions  allowed  in  that  region  (Figure  1).  Although 
there  is  a  change  in  sign  of  one  coefficient  at  A  =  A^,  there  is  no  corresponding 
change  in  the  class  type  of  solution  and  hence  A^  does  not  appear  in  the  figure.  A 
further  subdivision  of,  and  final  assignment  of  classes  in  the  regions  of  the  A,  B 
plane  is  accomplished  in  Section  4. 


Figure  1.  Division  of  (A,  B)  Plane  Into  Separate  Regions  Characterized  by 
Numbers  and  Classes  of  Solutions  in  Each  Region 
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We  now  proceed  to  examine  the  solutions  of  Eq.  (7)  so  that  classes  of  solutions 

may  be  assigned  to  the  different  regions. 

We  intend  to  examine  the  values  of  C  as  functions  of  A  for  fixed  <D  .  In  this, 

o 

we  define  Cj,  C ^  so  that  outside  of  values  for  C  at  the  singular  points 

At  =  (1  -  B)2  (B  +  1/2),  A2  =  (1  -  B)2  (B/2+  1)  we  will  have  C  *  C2  2  Cg.  The 
labelling  of  the  roots  in  this  order  also  identifies  them  permanently  by  branch.  If 
we  put  C=C'  -  p/3  then  Eq.  (18)  becomes 


(c')3  +  a  c'  +  b  =  0  . 

(20) 

«.  (-  b/2  +  V7T)1/3 

0  =  b/2-  Vr“  y3 

solutions  are  labelled  so  that 

c{  =  «  +  0 

(21a) 

C2=  - 2^  ^~""3 

(2lb) 

(21c) 

3.1  Asymptotic  limits  for  Small  A 

For  small  A  the  roots  of  Eq.  (7)  are  given  by  relatively  simple  expressions. 

If  the  equation  is  modified  so  that  terms  in  A  only  as  high  as  the  second  power  are 
retained,  we  have 


2C3 - <B2 


1>2  C2 


16/27  A(A(1+  B2)  -  <B2  ■ 


1)  (B3  -D)  C  + 


( 8 / 27 )2  A2 


(B3  -l)2  =  0 


(22) 


If  we  assume  a  solution  0|C]  >>  OjA]  where  Oj  ]  =  order  of  we  may  drop  the  term 
2 

in  A  ,  and  the  resulting  quadratic  has  as  one  solution 


C,  =  |  (B2  -  l)2 


16/27 


(23) 
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The  other  root  of  the  quadratic  leads  to  0[C]  =  0[A]  and  is  hence  contrary  to  our 

2 

assumption.  If  Eq.  (6)  is  expanded  out  to  terms  of  order  A  we  find 
V  2 C  =  ±  (B2-  1).  Thus  the  single  solution  of  Eq.  (22),  Cj  =  1/2(B2-  l)2  contains 
the  two  solutions  of  Eq.  (6)  corresponding  to  the  ±  sign. 

This  does  not  however  exhaust  the  solutions  of  Eq.  (22).  If  we  now  assume 

3 

that  there  is  a  solution  of  the  form  0(C]  =  0(A]  the  term  in  C  is  of  higher  order 
in  powers  of  A  than  the  other  three  terms  and  solution  of  the  quadratic  resulting 
from  neglecting  the  term  in  C3  yields,  for  B^  1 


C3  =  8/27 


A  =  8/27 


B  +  B  +  1 

irn 


(24) 


This  is  a  double  root  of  Eq.  (21)  and  hence  the  existing  three  roots  have  been  found. 

2  '1 

For  B=  1,  the  approximations  A  f(B)  «  (B  -  1),  (B  -  1)  with  f(B)#  0  are  no 
longer  valid  and  we  must  return  to  Eq.  (7).  Setting  B=  1  and  then  retaining  only 
the  lowest  power  of  A  in  the  coefficients  we  have 

C3  -  16/27  A2  C  +  2  ( 8 /  27 ) 2  A3  =  0  . 


The  solutions  are 

=  Cg  =  4/9  A  (double  root)  (25a) 

C3  =  -  8/9  A  .  (25b) 

If  Eq.  (24)  is  put  into  Eq.  (6)  and  the  roots  expanded  to  obtain  the  lowest  non¬ 
vanishing  order  in  A,  we  find  that  C2  is  not  a  solution  of  Eq.  (6).  Although  it  is  a 
solution  of  Eq.  (11a)  it  is  also  not  a  solution  of  Eq.  (10a). 

Thus,  of  the  three  solutions  of  Eq.  (22),  only  C(  is  a  solution  of  the  original 
problem;  and,  this  contains  the  two  solutions  corresponding  to  the  ±  sign. 

3.2  Other  special  Values  of  (’. 

For  A=Ap=(B+  I)3  (and  lt=  0),  Eq.  (7),  when  solved  by  use  of  standard  cubic 
formulae  yields  expressions  which  may  be  simplified  to  obtain 

Cj  =  1/6  (B+  l)2  (3  B24  14/3  B+ 3)  <26a> 

C2  =  C3  =  -  8/9  B(l  +  B)2  .  (26b) 
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Similarly,  for  A  =  A^~  (1  -  B)2  (B/2  +  1) 


Cl  =  1  / 18  (B  -  l)2(B2+2B+9)  (27a) 

C2  =  C3  =  2/9(B- 1)2B(B+ 2)  (27b) 

and  for  A  =  A1  =  (1-  B)2(B  +  1/2) 

Cj  =  1/18(B- 1)2(9B2+ 2B+ 1)  (28a) 

C2  =  C3  =  2/9(B- l)2  (2B+ 1)  .  (28b) 


FbrA  =  A^j  =  (B+  2)2  (B  -  1)  we  know  one  root  from  the  definition  of  A^;  namely 

=  2C  +  16/9  A.. 

^/y=0  M 

and 

C  =  Cg  =  -  8/9  (B+ 2)2  (B-  1)  .  (29a) 

From  this  we  obtain  an  equation  quadratic  in  C,  and  the  remaining  two  roots 
Cj  2=  1/4  |(B2-1>2+  16/9  (B+  2)2(B-  1)  ±  (B2-I) 

£(B2-  l)2  +  32/3  (B+  2)2  (B  -  l)j  W2|  (29b) 

where  the  ±  signs  are  associated  with  the  subscripts  1,  2  of  C  respectively. 
Similarly  for  A  =  A^  =  (2  B+  l)2  ( 1  -  B)  one  root  is 

Cg  =  -  8/9  A  B  =  -  8/9  B( 2  B+  l)2  (1  -  B)  (30a) 

and  the  other  two  roots  are  found  to  be 
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c 


1. 


2 


=  1/4 


|  (B2-  l)2  +  16/9  AB  ±  £  [(B2  -  1>2  +  26/9  Ab] 
-  64/81  AB(B  -  1)  (-25  B3  -  25  B2  +  5  B+  9) 


where  again  the  ±  signs  go  respectively  with  the  subscripts  1,  2. 

3.3  The  Zeros  of  C 

The  zeros  of  C  are  obtained  when 

4  A  -  (B3  -  1  -  A)2  =  0 


(30b) 


or 


A  =  A^  =  (B3/2  ±  l)2  .  (31) 

For  any  B  there  are  thus  two  values  of  A  for  which  C=  0.  It  must  be  deter¬ 
mined  however  which  branches.  Cj,  C'2  or  C 3  these  zeros  belong  to.  This,  as 
well  as  the  further  properties  of  C  require  more  complete  solutions  for  t:  as 
functions  of  A.  We  now  turn  to  these  solutions. 

3.4  Systematic  Presentation  of  Solutions  for  0 

We  continue  our  examination  of  solutions  of  Eq.  (7)  bv  graphing  solutions  of 

C  vs  A  for  seven  values  of  ip  :  -0.  9;  -0.  5;  0.  0;  1.  0;  3.  0;  30;  and  1000.  These  are 

o 

shown  in  Figures  2  through  8.  The  first  five  values  are  chosen  to  fall  in  the 
separate  regions  illustrated  in  Figure  1;  the  last  two  to  show  behavior  at  large  <PQ. 
In  addition.  Figures  9,  10  and  11  show  the  behavior  of  the  roots  for  small  A,  for 
<P  =  -0.  5;  1.  0;  and  3.  0. 

°  3/2  2 

We  note  first  that  the  zero  at  A  =  (B  -  1)  occurs  in  the  C.,  branch,  while 

3/2  2  J 

that  at  A  =  (B  +1)  occurs  in  the  C2  branch.  Further,  again  aside  from 

singularities  at  A  =  A^  and  A=  A2  the  intent  to  label  the  branches  so  that 

Cj  *■  C2  ?  Cg  is  accomplished.  We  next  turn  to  a  discussion  of  these  singularities. 
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Figure  4.  Grap.i  of  Integration  Constant  C  vs  Non-dimensionalized  Current 
Density  A,  for  Non-dimensionalized  Voltage  <PQ  =  0.  0 


°o*l  0 


0  2  4  6  0  10  a  14  16 


A 

Figure  5.  Graph  of  Integration  Constant  C  vs  Non-dimensionalized  Current 
Density  A,  for  Non-dimensionalized  Voltage  <PQ  =  1.0 


Figure  7.  Graph  of  Integration  Constant  C  vs  Non-dimensionalized  Current 
Density  A,  for  Non-dimensionalized  Voltage  «5  =  30.0 


\ 


A 


Figure  8.  Grapn  of  Integration  Constant  C  vs  Non-dimensionali2ed  Current 
Density  A,  for  Non-dimensionalized  Voltage  <Pq  =  1000 


%*-05  (Oetoil) 


Figure  9.  Graph  of  Integration  Constant  C  vs  Non-dimensionalized  Current 
Density  A,  for  Non-dimensionalized  Voltage  =  -0.  5.  Details  for  small  A 
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3.5  Singularities  in  C 

The  singularities  of  C  occur  at  the  branch  points  given  by  R(A)=  0.  We  con¬ 
fine  ourselves  to  some  general  comments  pointing  out  several  of  the  more  obvious 
features  of  these  singularities. 

We  note  to  start  with  that  as  far  as  first  derivatives  are  concerned  a  visual 
inspection  of  Figures  2  through  11  indicates  that  for  B  <  1  and  Aj  <  Ag,  Cj  and  C2 
have  discontinuities  at  A  =  Aj,  while  C'2  and  C3  have  discontinuities  at  A  =  Ag.  For 
B  >  1,  and  A^  >  the  roles  of  Aj  and  are  interchanged. 

For  the  behavior  of  the  functions  themselves  Eqs.  (27)  and  (28)  yield  the 
following: 


B 

> 

n 

> 

A=  A2 

B  <  1 

c  <  c  = 
C1  C2 

C3 

C1  >C2  = 

C3 

B  >  1 

C1>C2  = 

C3 

C1  <C2  = 

C3 

B  =  1 

C1  =  C2  = 

C3 

ci  =  c2  = 

C3 

At  A  =  the  smaller  of  A^  and  A2  the  evidence  of  a  discontinuity  is  clear.  The  fig¬ 
ures  show  at  this  point  that  Cj  =  C'2  >  C^:  Eqs.  (21)  on  the  other  hand  have  when 
K=  0,  in  general,  C'2=  C'3  4  Cj.  There  thus  exists,  for  each  fixed  value  of  B  a 
singularity  of  C  at  A  =  the  smaller  of  A^  A2-  This  singularity  has  the  effect  of 
exchanging  the  neighboring  values  of  Cj  and  C'3.  Numerical  solutions  are  con¬ 
sistent  with  this  supposition.  Thus  we  have  for  example  for  B  =  2  the  following: 


A 

C1 

C2 

C3 

1.  99997 

1.  77782 

1.  77775 

0.  944434 

1.  99998 

1.  77780 

1.77776 

0.  944437 

1.  99999 

1.  77779 

1. 77777 

0.  944441 

2.  00000 

0.  94444 

1. 77778 

1.  777778 

2. 00001 

1.  77779 

1.77776 

0.  944448 

2.  00002 

1. 77780 

1.  77775 

0.  944452 

2.  00003 

1. 77781 

1. 77774 

0.  944455 

At  A=  the  larger  of  A  ,  A0  there  is  no  obvious  requirement  that  a  singularity 
exist,  since  graphical  values  in  the  figures  and  the  results  of  Eqs.  (21)  both  lead 
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3.6  Assignment  of  (lxstu1*  to  Solutions 

Having  obtained  the  solutions  of  Eq.  (7)  it  now  remains  to  determine  which  of 
these  satisfy  Eq.  (6)  (class  I)  and  which  satisfy  Eqs.  (10a)  and  (lib)  (class  II). 

This  is  accomplished  by  computer  by  substituting  values  of  C  into  the  appropriate 
equations.  If  this  is  done,  changes  in  numbers  and  classes  are  found  to  occur  only 
for  A=  A^j,  A^  .  and  Aj  .  This  information  is  presented  in  Table  1.  ihe 
assignment  of  solutions  in  class  III  are  discussed  in  the  following  section. 


Table  1.  Numbers  and  Classes  of  Solutions  in  the  Various  Regions 
of  Parameter  Space 


Region 

Solutions 

Class  I 

Class  11 

Class  111 

1 

1 

1 

0 

0 

2 

1 

1 

0 

0 

3 

1 

0 

1 

0 

4 

3 

0 

2 

1 

5 

3 

1 

1 

1 

6 

3 

1 

1 

1 

7 

1 

0 

0 

1 

8 

3 

1 

0 

2 

9 

3 

0 

1 

2 

3.7  Uniqueness 

For  A  <  At  (classes  I  and  II)  and  A  >  AD  (class  Ill)  the  solutions  for  C  are 

i-4  X  ** 

unique  functions  of  <0Q  and  A.  For  A^  <  A  <  A^on  the  other  hand  the  solutions  are 
many  valued.  Which  solution  among  the  multiple  solutions  is  actually  realized,  may 
be  determined  by;  (1)  the  manner  in  which  the  operating  conditions  are  approached, 
and  (2)  questions  of  stability  of  operation.  We  do  not  investigate  either  of  these 
questions  here. 

3.8  Determination  of  Values  of  C 

Value(s)  of  C  for  any  <Dq,  A  may  be  found  by  use  of  Table  2.  ihe  table  is  con¬ 
structed  by  examination  of  solutions  corresponding  to  the  regions  shown  in  Fig¬ 
ure  1  and  by  observing  that 
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A 

A 

A 


l"  Ar  A2 

-i<po 

M  >  Aj,  A2 

9° 

VI 

H 

>  A  A 

m  '  r  a2 

0  ±  <0 

o 

<  00 
£  0 
<  B 


Table  2.  Assignment  of  C  to  Branches  Corresponding 
to  Ranges  of  A  and 


-15<P  S  0 

O  L 

and 

%-«3o<C° 

(AI~  Al*  A2> 

Range  of  A 

Branch  of  C 

Class 

A  £  Aj 

C1 

1 

A1S  AiA2 

C  2 

I 

a2  5  A  - Al 

C3 

I 

a+l  £  A  £  A^ 

^2^3 

II/l 

amsasar 

C2'C3 

II 

°<AM2  A1 

.  A9)  and  0S«  < 

*  o  U 

(Am >  Ar  A2) 

Range  of  A 

Branch  of  C 

Class 

A  £  Aj 

C1 

I 

A1SASA2 

C2 

I 

a2  ~  A  ~  Am 

C3 

I 

am2A5A1 

C3 

11 

A  L,  ~  A  ”  A  R 

C2'  C3 

II 

Ati  = 
A*r 

At,  for  (p  >  0 
M  o 

A. ,  for  ip  <  0 
M  o 

THE  SPACE  CHARGE  LIMITED  (CLASS  111)  CASE 
We  introduce  the  following  additional  terms  (Kigure  12); 

,j  =  Fractionof  emitted  current  which  is  transmitted  past  the  potential  minimum, 
y  j  =  Value  of  y  at  the  potential  minimum. 
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<p=  -1  +  (A  j)2/3  (y  -  y  j)4/3 


respectively.  The  boundary  conditions  at  y  =  0  and  y  =  1  then  yield 


1  =  A(2-j)yJ 


B3  =  Ajil-yj)2  . 


Eliminating  j, 

2  A  =  l/y2  +  B 


(35) 

(36) 

(37) 

(38) 


d-yj) 


and  clearing  of  fractions. 


F  =  2  Ay4  -  4  Ay3  +  (2  A  -  1  -  B3)  y2  +  2y  -  1=0 


(39) 


an  equation  of  fourth  order  in  y^  The  role  of  yt  here  replaces  that  played  by  the 
integration  constant  in  the  other  two  cases;  once  it  is  known  (a  function  of  A  and  B) 
the  potential  distribution  is  fixed. 

Since  F(0)  =  -  1  and  F(-o0)  =  <*>,  Eq.  (39)  has  at  least  one  real  root  between  0  and 
3 

Since  B  C 1 )  =  -B  and  F (°°  ) = » ,  there  is  at  least  one  real  root  between  +1  and 
+  °°.  F rom  Eq.  (38)  it  is  clear  that  the  other  two  roots,  when  real,  lie  between 
yt  =  0  and  y  j  =  1 . 

Consider  Eq.  (38)  for  a  fixed  value  of  B.  For  y  j  =  1  or  y  ^  =  0,  A  .  For 
some  0  <y^  <1,  A  has  a  (unique)  minimum.  A  must  always  be  larger  than  this 
minimum,  and  finding  this  minimum  supplies  a  necessary  condition  in  defining  the 

domain  of  A,  B  space  over  which  there  will  be  allowable  solutions  for  v..  We  have. 

J  1 


dA 


-2  .  2  BJ 


'  3yT  ~  *  - - 3 

1  yi  d-yj) 


/>  u3  3 

(1  -  yj)  =  B  yj 


=  0 


(40) 


of  physical  significance  exists  for  A  >  A  j 
by  graphical  display  in  Figures  13a  and 


.  .  These  considerations 
13b. 


Figure  13a.  Geometric  Illustration  of  Limits  and 
Numbers  of  Class  III  Solutions  for  B  2  1 


Figure  13b.  Geometric  Illustration  of  Limits  and 
Numbers  of  Class  III  Solutions  for  B  §  1 


We  note  that  A^j^  intersects  Aj^  at  B  =  1  I*  >  B^  and  intersects  A^  at  B=  1. 
Figure  14,  a  superposition  of  A^j^  onto  Figure  1  completes  the  classification  of 
solutions,  dividing  the  A,  B  plane  into  nine  distinct  regions.  The  region  labelled 
2  in  Figure  1  is  split  into  two  regions,  2  and  8,  in  Figure  14;  and  the  region 
labelled  3  in  Figure  1  is  likewise  split  into  regions  3  and  9  in  Figure  14.  The 
numbers  of  class  III  solutions  in  each  region  has  already  been  given  in  Table  1. 
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Figure  14.  Division  of  <A,  B)  Plane  Into  Separate  Regions  Characterized  by 
Numbers  and  Classes  of  Solutions  in  Each  Region.  Detail  for  B  5  1  shows 
modification  of  Figure  1  to  include  class  III  solutions 


4.1  Special  Values  of  vj 

For  several  particular  values  of  B,  the  solutions  of  Eq.  (33)  are  of  simple 
form. 

For  B=  0,  the  four  roots  are 

y,  =±  ==  ,  1  (the  last  a  double  root)  . 

1  }/2A 

For  U=  1,  the  four  roots  are 


=  ± 


A  +  2  ± 


2  V2  A  +  1 

TA - 


1/2 


+  2 


The  two  roots  0  "  *  1  are  those  for  which  the  ±  sign  inside  the  radical  is  taken 
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as 


For  A  =  (B 


7 


3/2 


+  1)  one  root  is 


yl  = 


Va 


The  residual  cubic  is 

2yfAy3l  +  (2  -  4  VA)  y2  +  (Va  -  2)  y  {  +  1  =  0 

and  the  remaining  roots  do  not  seem  to  have  simple  forms. 
For  A  =  1  - —  ,  the  four  roots  are 

=  Y^-g  (double  root) 

_  B  ±  Vb2  +  B  +  1 
yl  "  B  +  1 


Figure  15  shows  typical  values  of  y^  as  a  function  of  A  with  B  as  parameter. 
Solid  lines  are  solutions  of  the  original  boundary  value  problem  with  j  “=  1.  Dashed 

values  continue  solutions  for  all  values  Oiy.  s  1.  Dotted  lines  (1/dA  for  B  5  1 

—  -1/3  ^  " 

and  1/VA  and  (2  A)  for  B  S  1)  separate  regions  according  to  number  of  class  111 

solutions  existing  in  each  region. 


Figure  15.  Values  of  0  s  y,  <  1  From  Eq.  (33)  as  Functions  of  A  With  B  as 
Parameter.  Solid  lines;  values  of  yj  for  which  j  -  1;  dashed  lines;  values  of 
y.  for  which  j  -  1;  dotted  lines;  separate  regions  according  to  number  of 
class  III  solutions  existing  in  each  region 
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FORMS  OF  THE  POTENTIAL 


The  primary  purposes  of  this  report  are  to  investigate  structure,  to  classify 
solutions,  and  to  organize  methods  of  obtaining  the  latter.  We  do  not,  therefore, 
present  here  large  numbers  of  graphs  of  potential  vs  distance;  these  are  given  in 
works  previously  referred  to.  We  do,  however,  for  illustrative  purposes,  give 
several  examples  of  such  graphs  and  to  this  end  we  choose  parameters  lying  in 
regions  4,  5,  8  and  9,  in  all  of  which  there  are  three  solutions  corresponding  to 
each  pair  of  parameters  (A,  B).  Figures  16,  17,  18,  and  19  each  show  the  three 
potential  curves  corresponding  to  the  values  of  parameters  shown  in  Table  3. 


Table  3.  Parameters  For  Graphs  of  Potential.  Numerals  in  Parenthesis 
Give  Class 


Figure  Number 

16 

17 

18 

19 

Region 

4 

5 

8 

9 

B 

1.  0 

2.  0 

0.  1 

0.  6 

o 

o 

0.  0 

3.  0 

-0.  99 

-0.  64 

A 

6.  0 

15 

1.  0 

2.  1 

C1 

... 

... 

--- 

... 

C2 

-0.  9261(11) 

-0.  1651(11) 

... 

... 

C3 

-4.  086  (II) 

-13.31  (I) 

-0.  02545(1) 

-1.  118(11) 

*1 

yl 

0.  3187 

0.  2525 

0.  7092 

0.  9672 

0.  5794 

0.  6679 

j 

.i 

0,  3591 

0.  9546 

0.  0118 

0.  9311 

0.  5813 

0.  9324 
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Appendix  A 

Listing  o'  Special  Value  of  A 


M 


M 


kR 


A 

A 


_  (b2,1>1B3-D(B  +U 

e  (2  B+  l)2  ^ 

=  (B+l>3 


-1 


=  i(B+i> 

MIN  2  o 

=  U-B)2  (B/2+D 


(Aj.  A2*  AR 


are 


roots  of  the 


discriminant  of  Eq. 


(7).> 
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